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ABSTRACT 

We show that when the Hornberger-Sipe calculation of collisional 
decoherence is carried out with the squared delta function a delta of energy 
instead of a delta of the absolute value of momentum, following a method 
introduced by Diosi, the corrected formula for the decoherence rate is simply 
obtained. The results of Hornberger and Sipe and of Diosi are shown to be 
in agreement. As an independent cross-check, we calculate the mean squared 
coordinate diffusion of a hard sphere implied by the corrected decoherence 
master equation, and show that it agrees precisely with the same quantity 
as calculated by a classical Brownian motion analysis. 
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1. Introduction 

The calculation of collisional decoherence was initiated by Joos and Zeh [1], with 
generalizations of their result and corrections to the overall normalization given in papers of 
Gallis and Fleming [2] , Dodd and Halliwell [3] , and Hornberger and Sipe [4] . An even more 
general calculation of collisional decoherence was also given by Diosi [5], but this was not 
known to Hornberger and Sipe, while noted by Dodd and Halliwell, as well as in the master 
equation papers of Altenmiiller, Miiller, and Schenzle [6] and of Vacchini [7] . A difficulty 
encountered in refs [1], [2], and [4] is the appearance of a squared delta function of the 
absolute value of momentum in the calculation of the decoherence rate. To circumvent this, 
Hornberger and Sipe carried out a careful wave packet analysis, which gives an answer smaller 
by a factor of 2n than that given by Gallis and Fleming, and this result is in agreement with 
experiment [8]. 

In the course of an alternative derivation of the corrected result, following the "tra- 
ditional approach" of refs [1] and [2], Hornberger and Sipe introduce a rule in which the 
squared delta function of absolute value of momentum is evaluated in terms of an inverse 
scattering cross section, which drops out later in their calculation. This mixing of kinematic 
quantities (such as a squared delta function) and dynamical ones (such as a cross section) 
is unconventional, and Hornberger and Sipe describe this part of their calculation as specu- 
lative. In Sec. 2 we show that an entirely conventional completion of the Hornberger-Sipe 
calculation is possible, if one follows the method used in the earlier and more general master 
equation derivation given by Diosi, and also used in the derivation of Dodd and Halliwell. 
These authors retain the delta function of energy that appears as the T-matrix coefficient, 
rather than converting this delta function to a delta function of the absolute value of mo- 
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mentum. This makes a difference when squaring the delta function. For a squared delta 
function of energy, one can use the standard rule of evaluating 5(0) in terms of the elapsed 
conjugate time variable, as is done in the usual textbook "golden rule" calculation. For a 
delta function of absolute value of momentum, it is not so clear what to use as the corre- 
sponding conjugate variable when taking the square, and this appears to be the root of the 
difficulties in the earlier calculations of refs [1], [2], and [4]. 

In Sec. 3 we compare the Hornberger-Sipe and Diosi results and show that they are 
the same; hence Diosi's 1995 calculation appears to be the first giving the correct result for 
the collisional decoherence rate. In Sec. 4 we give an independent check of the corrected 
expression for collisional decoherence, by using the corresponding master equation to evaluate 
the scattering-induced translational Brownian diffusion of a hard sphere in the geometric 
scattering limit. Planck's constant cancels out in this calculation, and so the result obtained 
this way can be directly compared with the classical Brownian diffusion of a hard sphere, 
and the two calculations are in precise agreement. 

2. Calculation using a squared energy delta function 

To keep this section concise, we will use the notation of Hornberger and Sipe, and 
give just a brief summary of their calculation up to the point where our treatment begins 
to differ from theirs. We consider a Brownian particle in a bath of N scattering particles 
of mass m, contained in normalization volume Q. In the dilute case, the scatterings of the 
bath particles from each other can be neglected, and their scatterings from the Brownian 
particle are independent of one another. Then the effect of the N bath particles is obtained 
by considering the effect of a single bath particle, and multiplying by iV at the end of the 
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calculation. For a single scattering, the effect of the collision is to change the density matrix 
p (Ri, R2) to 

p(Ri,R 2 ) =r ? (R 1 ,R 2 ) P o(Ri,R 2 ) , (1) 
with the factor ?y(Ri, R 2 ) given by 

?7(Ri,R 2 ) =tr bath {e^ p - R2/?i l S t e ip - {R2 - Rl)/n l Soe ip - Rl/ V ath } • (2) 

Here So is the scattering matrix, and p bath corresponds to an ensemble momentum space 
weighting 

j dpnip) , (3a) 
where for a thermal ensemble with (3 = (kT)^ 1 one has 



2irm / 

Substituting So — 1 + i%, evaluating the bath trace in a momentum basis, inserting a 
complete set of intermediate states, changing from box to continuum normalization, and 
using the unitarity relation i(T — Tq) = —TqT , Hornberger and Sipe show that Eq. (2) 
takes the form 



^(Ri,R 2 ) = J d W (p) [l - / dp'(l -e 4(p - p ' HRl - R2)/?i )|(p'|To|p) 



(4a) 



This equation, which is Eq. (51) of Hornberger and Sipe, and without the weighting over 
yu(p) also corresponds to the first line of Eq. (2.11) of Gallis and Fleming, will be the starting 
point for our analysis. 

Denoting the elapsed time in the scattering process by T, so that Eqs. (1) and (4a) 
implicitly refer to time T, and recalling that 1 = J dp/i(p), we can rewrite Eqs. (1) and (4a) 
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as 

p(Ri, R 2 ; T) - p{R u R 2 ; 0) = (^(Ri, R 2 ; T) - l)p (Ri, R 2 ) 
= - Po (Ri,R 2 ) / dMp)^f- 1 rfp'(l-e^'H R - R ^)|(p'|T |p)| 2 . m 
Our next task will be to evaluate the squared matrix element appearing in the integrand of 
Eq. (4b). 

The general T matrix element (cfcl^olqi) can as usual be expressed in terms of the 
scattering amplitude /(q 2 ,qi) and an energy-conserving delta function, 

(qal^oloi) = ^f~ 6 ( E2 ~ s i)/(q2, qi) 

\ m (5) 
= 2^ fe - ?l)/(q2 ' qi) • 
Instead of using the second line of Eq. (5) to form the square of the % matrix element, 

we will use the first line, using the second line only after the delta function of zero energy 

argument has been evaluated. Thus, we have 

l(q2|To|q 1 )| 2 = ^^5 2 ( J E 2 - J E 1 )|/(q 2 ,q 1 )| 2 . (6) 

Using the Fourier representation for the energy delta function, 

1 f°° 

5(E 2 -E 1 ) = —j dtexp[t(E 2 - E 1 )t/h] , (7a) 



we find 



5 2 (E 2 -E 1 )=5(E 2 -E 1 )5(0) 



dt 



m T Rl \ 



(76) 



q 2 2nh 

where in the final line we have converted the energy delta function to a delta function of 
the absolute value of the three-momentum. Here T is the elapsed time interval, which we 



assume to be longer than the time for a single scattering, but still short compared to the 
characteristic decoherence time of the Brownian particle in the iV-particle bath. 

Replacing q 2 ,qi by p', p respectively, substituting Eqs. (7b) and (5) into Eq. (4b), 
and writing dp' = dn(p') 2 dp', with dn a solid angle differential, we can immediately integrate 
out the delta function of the absolute value of momentum. Using 

(2nh) 3 1 2 m T T p 

TP —7T^ = o— ' ( 8a ) 



Vl (2irhm) 2 p 2nh f2 m 



we thus get 



p(Ri, R 2 ; T) - p(Rx, R 2 ; 0) = - p (Ri, R 2 

T 



rfp/i(p)£ J dn(l - e i(p - pftHRl - R2)/?i )|/(pn,p)| 2 . 

(86) 



Multiplying by iV to take account of the fact that each of the N bath particles makes 
a contribution equal to Eq. (8b), denoting the bath density N/fl by n, dividing by T, and 
finally passing to the limit of small T, we get the result 

9p(Ri , R 2 ; t) = _ _ R 
at 

with 

F(Ri - R 2 ) = n j dp/i(p)£ J dh{\ - e i(p - pft) - (Rl - R ^ )/ ' i )|/(pn, p)| 2 . (96) 

This is the form of the final result for the decoherence-induced master equation given in 
Eq. (55), and in the unnumbered immediately preceding equation, of Hornberger and Sipe 
[4]. By defining u(p) by 

Kp)dp= <M^i , (10 .) 
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(10c) 



with a second solid angle differential, so that J °° dpu(p) = 1, Hornberger and Sipe also 
rewrite Eq. (9b) in the equivalent form (after a relabeling of the integration variables) 

F(R) = n r dqv(q)- f (1 - e iq ^-™W n )\f{qh 2 , qn^ 2 . (106) 

Jo m J An 

In the limit of large R, the exponential term in Eqs. (9b) and (10b) averages to zero 

provided that n 2 ^ "hi, while in the forward direction h 2 = "hi the integrand in Eqs. (9b) 

and (10b) vanishes for all R. Hence one has 

F(R->oo)=n/ dqv(q)— [ < ^^\f(qn 2 ,qh l 
Jo rn J h2 ^ hl 4vr 

=n / dqv(q)—<j(q) , 
Jo m 

with a(q) the total cross section (excluding a possible delta function contribution to the for- 
ward diffraction peak). In other words, the large R asymptote of F(R) is the thermal ensem- 
ble average (nva) av of the non-forward scattering rate n(q/m)a(q). Correspondingly, from 
Eqs. (4b) and (8b), as modified by multiplication by the factor N, we see that ^(Ri, R 2 ; T) 
= 1 — TF(Ri — R 2 ) has the large R asymptote 77(00; T) = 1 — T{nva)p^. Thus 77(00; T) 
vanishes for T equal to the inverse of the averaged non-forward scattering rate, a result 
reminiscent of, but not identical to, the condition 77(00) = imposed by Hornberger and 
Sipe on the single collision decoherence function in their Eq. (45), on which they base their 
method for evaluating the square of a delta function of the absolute value of momentum. 
(We emphasize, however, that in the calculation leading to Eq. (9b) we have not fixed T by 
imposing such a condition.) 

3. Comparison of the Hornberger-Sipe and Diosi Results 

The calculation of Diosi includes effects of recoil of the Brownian particle; we show 
in this section that in the limit of an infinitely heavy Brownian particle, the results of 



Hornberger-Sipe and of Diosi are in agreement. When recoil is neglected, Diosi's result is 
his Eq. (19), which reads 

f t = n J dEd^dttfe^p-p^k,) (v^pVl^ - \{Vi } ^ v p}j . (11a) 

Substituting into Eq. (11a) Diosi's Eq. (20) (with the Brownian particle recoil term dropped), 

Vkjki = exp(-ik /i • q) = exp(-i(k f - k;) • q) , (116) 

taking the matrix element of Eq. (11a) between q eigenstates (Ri| and IR2), and writing 

p(Ri,R2;f) = (Ri|p(*)|R2) , (11c) 

we get 

^ (Rl rf f 2;t) = -no f dEdQ t dn f e d -^lp £ (k t ) (1 - e^-M-CRi-*.)) p(Rl , R 2; t) . 

(12) 

Taking account of the fact that dE = dk 2 /(2m) = (k/m)dk, together with da(9 : E)/dflf = 
\f\ 2 , where / is the scattering amplitude, along with p £ (kj) = /i(kj) and some obvious 
relabelings of variables, one sees that Eq. (12) is identical to Eqs. (9a) and (9b) that follow 
from the analysis of Hornberger and Sipe. 

Diosi also gives the expansion of Eq. (12) to leading order in R = Ri — R 2 . Working 
now in the other direction, from Eq. (10b), the simplest way to find the leading order R 
dependence is to note that F(H) is a rotationally invariant function of R. Hence it suffices 
to evaluate the average over the direction of R. Expanding out the exponential in Eq. (10b), 
we have 

x _ ^(m-myn/n = x _ _ ^ . R / n + htf /K 2 )[(n x - h 2 ) ■ R] 2 + ... . (13a) 
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The average of [(hi — h 2 ) • R] 2 over the direction of R is 



i i [i/t / / 

3 z/ 1 3 

(136) 

=|i2 2 (l-ni.n 2 ) = |i2 2 siii 2 (0/2) , 
with the angle (the scattering angle) between hi and h 2 . Substituting Eqs. (13a,b) into 

Eq. (10b), we get 

F(R) = R 2 A , (14a) 



with 



. 2 n f°° , , . q 2 f dh\dfi2 . 2//1/ „ x , ,, * 

A= oT2 / dqv(q)—q / — sm (9/2) \f(qn 2 , qn x 

3 h z J m J An 

2 n f q f 

= 3^2 / dqfi(q)—q 2 / dh 2 sm 2 (9/2)\f(qh 2 , qh 



(146) 



With the shifts in notation noted above, this equation for A is identical to Eq. (22) of Diosi, 
which gives what he terms the diffusion parameter D pp . 

4. Comparison of decoherence-based and classical 
calculations of Brownian translational diffusion 

Let us now apply the result of Eqs. (14a,b) to the quantum Brownian motion of a 
Brownian particle of mass M and radius a, with a large enough so that pa » 1 for important 
bath particle momenta. In this case the quantum scattering differential cross section consists 
of two parts [9]: an isotropic part, with an integrated cross section contribution of na 2 , 
and a forward diffraction peak, again contributing cross section ira 2 . Since the integrand of 
Eq. (14b) vanishes for forward scattering, the forward diffraction peak makes no contribution. 
Thus we can evaluate the integrals by taking the scattering amplitude to be a constant, 
f(qh 2 ,qh 1 ) = F, with 

r ds\F\ 2 = 4n\F\ 2 = na 2 , (15a) 
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that is, with 

\F\ 2 = a 2 /A . (156) 

Substituting into Eq. (14b), we then find that 

A= nna^U (16) 

Here we have defined 

(q 2 v)AV= dqu(q)q 2 ^=A(m/n)H2kT)i , (17a) 
Jo m 

where we used the expression 

3 

u(q) = W (J- V e-WlW , (176) 
\27rm / 

which follows from the definitions of Eqs. (3b) and (10a), to evaluate the thermal average. 
Substituting Eq. (14a) into Eq. (9a), we then get 

dp(R^R 2 ; t) = _ A(Ri _ R2)2p(Ri? R2 , f) (lgo) 

which with Eq. (11c) is equivalent to the operator equation 

dp(t) 



-Aj2[Rj,[Rj,P(t)}} , (186) 



dt 

with Rj denoting the Cartesian components of R. Adding the kinetic energy term to the 
differential equation for p(t), we get the total evolution equation 

^ = -^[H km ,p(t)]-Aj2[R 3 AR 3 ,P(t)}} , (19a) 

3=1 

with the kinetic Hamiltonian for the Brownian particle given by 

H ^ = tw • ( 19i >) 

j'=i 
n 



Here Pj is the momentum operator corresponding to the coordinate operator Rj, so that 
[Rj,Pk] = iftSjk- 

From Eqs. (19a) and (19b), one can calculate the mean squared coordinate diffusion 
as a function of time, for a Brownian particle that starts at R = at t = with zero drift 
velocity. This calculation in the one-dimensional case is given in Adler [10], by constructing 
a generating function for the trace of pit) multiplied by an arbitrary polynomial constructed 
from Rj and Pk- For the mean square coordinate deviation, the result on converting to the 
present notation is (with no sum implied over j) 

2\fr 2 t 3 

(R])=tr P m*= 2 i^r ■ (20a) 

Substituting Eq. (16) for A and Eq. (17a) for the thermal average of q 2 v, we get finally (again 
with j unsummed) 

(R 2 .) = C(kT)lnm^a 2 t 3 /M 2 , (206) 

with 

C = ^(2n) l > . (20c) 

Note that the Planck constant H has dropped out of this result. Hence the formula of 
Eqs. (20b, c) is a classical result, and should be recoverable by a purely classical calculation. 

The formulas needed for a classical Brownian motion evaluation of (R 2 ) are sum- 
marized in a recent paper by Collett and Pearle [11]. Their Eq. (2.2) gives (again with j 
unsummed) 

<^> = ^ • < 21 «> 

with £ a viscosity factor, which for a sphere of radius a in a dilute bath is given by their 
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Eq. (2.5), 

o 

f = -na 2 (2irmkT)^ . (216) 
3 

Substituting Eq. (21b) into Eq. (21a) then gives a result identical to Eqs. (20b, c) above. 
Since the result of Eqs (20b, c) is directly proportional to the normalization constant in the 
collisional decoherence rate, this agreement gives added confirmation of the correctness of 
Eq (9b). 
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